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ABSTRACT 

We estimate the thin disk scale height of the Galactic population of L- & T-dwarfs based on star 
counts from 15 deep parallel fields from the Hubble Space Telescope. From these observations, we have 
identified 28 candidate L- & T-dwarfs based on their (i' — z') color and morphology. By comparing 
these star counts to a simple Galactic model, we estimate the scale height to be 350 ±50 pc that 
is consistent with the increase in vertical scale with decreasing stellar mass and is independent of 
reddening, color-magnitude limits, and other Galactic parameters. With this refined measure, we 
predict that less than 10 9 M Q of the Milky Way can be in the form L- & T-dwarfs, and confirm 
that high-latitude, z~6 galaxy surveys which use the i'-band dropout technique are 97-100% free of 
L- & T-dwarf interlopers. 

Subject headings: stars: low-mass, brown dwarfs — Galaxy: structure — Galaxy: stellar content 



1. INTRODUCTION 

The method of counting stars to infer the shape and 
size of the Galaxy is as old as astronomy itself. Earli- 
est efforts using this technique were famously flawed as 
they often relied on insufficient data or w holly incorrect 
assumptions (eg. lHerscheUll785t lKaptevdll922ft . Aided 
by ad vanced technology. iBahcall fc Soneiral fBfcS: 119801 
1981) demonstrated that the true power of star counts 
is realized when they are compared to simul ations of the 
funda mental equation of stellar statistics ijvon Seelieerl 
1898). The B&S method relies heavily on the assumed 
luminosity functions and density distributions, and has 
been a standard method for many subsequent studies. 

The Galaxy is traditionally characterized by having a 
Population I disk and Populatio n II spheroid. In a series 
of st udies of the exp onential disk. lGilmore & Reidl l|1983f) 
and iGilmorei l) 19841) established the need for a thick and 
thin disk whose scale heights are inversely pr oportional to 
the m asses of the studied stars (see Table 1 in lSieeel et alJ 
2002) . The standard description of the Galactic halo is a 
de Vaucouleur or power-law profile, while the Besancon 
flattened spheroid with c/a~ 0.5— 0.8 llBahcall fe Soneira 



1984 lRo"hm. Revle. fc Crezel2000t 171 arson & Humphreys 
2003; Ro bin et al.ll2003l) being the currently favored pa- 
rameters. Thorough discussions of star counts and their 
relevance to Galactic struc t ure are given in the Annual 
Review s bv iBahcalll (|198fiD iGTlmore. Wvse. fc Kuiikenl 
(1989J), and lMaiewsHl|1993|) . 

Many of the Galactic models and the majority of 
the literature examine on relatively luminous dwarf 
and/or giant stars and rarely address sub-stellar objects. 
The discover y of the first extra-solar, sub-stellar object, 
Glies e 229B ijNakaiima et al.lll995t iQppenheimer et ail 
Il995|) motivated the creation of the L and T spectral 
classes. With surface temperatures ranging from 750- 
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2200 K l|BnrgasserHl999lh the L- & T-dwarfs can con- 
taminate searches for z — 6, i' -band dropout objects 
ijYan. Windhorst, fc Coher]l2003j) by mimicking the ex- 
tremely red broad-band colors. This effect has remained 
largely unquantificd due to insufficient knowledge of the 
L- & T-dwarf IMF, Galactic distribution, and local num- 
ber density. P revious work on their IMF and local 
number density <|Reid et alJll999T: iChabrierl 120011 120021 
iLiu et aiT2 002) has suffered from limited statistics. With 
the deep imaging of 15 Hubble Space Telescope (HST) 
Advanced Camera for Surveys (ACS) parallel fields, this 
study increases the number of faint dwarfs by adding 
28 to the growing list. Little wo rk has be e n don e on 
the scale height of L- & T-dwarfs; ILiu et all (|2002f) and 
IPirzkal et alJ l|2005j) estimated 100-400 pc as based on a 
single object or a 3-4 objects in a single field. Therefore 
the primary goal of this work is to estimate the scale 
height of the L- & T-dwarf population by comparing 
the surface densities from the ACS parallel fields to the 
Galactic structure models. 

2. OBSERVATIONS 

The L- & T-dwarf candidates were selected from 15 
HST/ ACS parallel fields covering a broad range in Galac- 
tic latitudes and longitudes (see Tabled . All ACS fields 
have at least three independent exposures in F775W 
(SDSS-i') and F850LP (SDSS-z') with a total exposure 
time of 2-10 ks per bandpass . All fields are >90% com - 
plete at z'(AB)~26.0 mag (jYan fc Windhorst 2004a). 
We adopt the AB magnitude system iJOke fc GunrJ 
1983). 

After combining the individual ACS frames into fi- 
nal stacks using the PyRAF-based script multidriz- 
zle lIKoe kemoer et alJ I2002j) . the SExtractor package 
ijBertin fc Arnoutslll996jr was used in double-input mode 
to perform the matched-aperture photometry. The 
F850LP stack was used to define the optimal apertures 
for the flux measurements in both stacks. For source de- 
tection, we used a 5 x 5 Gaussian smoothing kernel with 
a FWHM of 2.0 pixels, which is approximately the same 
as the FWHM of the ACS point-spread function (PSF) 
on both image stacks. We used total magnitudes (corre- 
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Fig. 1. — Color-magnitude diagram for all point sources detected 
in our 15 ACS fields. The small dots represent all objects that 
met the stellar morphology classification, the solid stars are the 
28 L- & T-dwarf candidates, and the dashed lines represent the 
imposed color-magnitude limits. Many point sources associated 
with diffraction spikes, field edges, and spurious detections (objects 
smaller than the PSF, likely residual cosmic rays) brighter than 
z' = 26 mag were manually removed. The similar points fainter 
than z' = 26 mag were not removed from this figure. The sample 
of stellar candidates, becomes incomplete for z' > 25 mag. 
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Fig. 2. — The observed FWHM for all objects in Field 1 is plotted 
as a function of apparent z' magnitude. The small dots, asterisks, 
and filled stars represent all objects, the point sources used to define 
the "stellar-locus," and the objects selected as candidate L- & T- 
dwarfs. Clearly many objects were recognized as point sources, 
however only two met the (i'—z') color criterion. 



sponding to the MAG-AUTO option in SExtractor) for 
the photometry and adopted the zero p oints published in 
HST ACS Instrument Science Report (|De Marchi et all 
2004). 

As a cursory selection, all objects with (i'—z')>l.3 mag 
and z' <26.0 mag were considered. Figure ^ is the 
color-magnitude diagram for all point sources detected 
in all 15 fields. Dashed lines indicate the imposed color- 
magnitude limits and the large, filled stars represent the 
candidate L- & T-dwarfs. For the five fields for which the 
F475W (SDSS-(/) band was available, we required candi- 
dates to have (g' — i')>0 mag. Objects near the detector 
edges were not considered, yielding an effective area of 
~9 arcmin 2 per ACS field. Extended objects were elimi- 
nated from the analysis by using the FWHM parameter 
of SExtractor. In FigureEl we plot the FWHM as a func- 
tion of apparent magnitude for each object in Field 1 as 
small dots. The locus of points at FWHM~ 0'.'13 and 
z' < 26 mag (hereafter the "stellar locus"), represented 
as asterisks, are the unresolved objects with the mini- 
mum possible FWHM. In addition to the above color- 
magnitude criterion, we required all L- & T-dwarf candi- 
dates to lie within this locus. Sources of contamination 
to these criteria are discussed in § 14.11 

Typically Galactic structure studies examine star 
counts from on e or many shallow fields with large sur- 
veyed areas (eg. iSiegel et al.ll2002t lLarsen fc Humphreys! 
20031) Thus the data of the 2M ASS iBurgasserl 
T^lKrrkpatrick et alJll999l). DENIS iDelfossd 119991) . 
and/or SPSS llStrauss et alJll999l: l Tsvetanov et alJl2000l: 
lHawlev et alJ I2002D are natural choices to study the 
Galactic distribution of the L- & T-dwarfs. These sur- 
veys have significantly more detection area than our HST 
dataset and their L- & T-dwarfs are typically closer to 
the Sun than ~300 pc or 1 disk scale height. Therefore 
to avoid extrapolating the vertical scale height beyond 



TABLE 1 
T-Dwarf Number Counts 



Field 


I 11 


6 U 


Number 


No. 


(deg) 


(deg) 


(per 9 arcmin 2 ) 


1 


115.018 


+46.681 


2 


2 


164.056 


-75.750 


1 


3 


169.188 


-59.664 





1 


279.934 


-19.990 


3 


5 


165.876 


+36.396 





6 


280.782 


+68.293 





7 


113.107 


+28.548 


2 


8 


293.996 


-41.466 





9 


316.829 


-40.490 


1 


10 


105.103 


+7.075 


11 


11 


298.138 


-13.885 


6 


12 


92.666 


+46.378 


1 


13 


70.106 


+62.876 





14 


251.327 


-41.444 





15 


216.142 


+54.561 


1 



this distance, we chiefly analyzed the HST dataset where 
all sample stars are likely more distant than 1 e-folding 
length. 

3. THE SIMPLE GALACTIC MODEL 

The Galactic structure models were made by distribut- 
ing 10 10 points according to an exponential disk that 
was motivated by the light profiles of ed ge-on galaxies 
ijde Griis. Peletier. fc van der Kruitj fl996l) : n(r,9,z) — 

n exp ( i ^f r ) exp (j^A , where r z =2100 pc is the 

radial scale length found by IPorcel et all (l998), 
n =0.12 pc~ 3 is the loca l spac e density L- & T- 



dwarfs take n from iChab rierl <l20 02f) . a nd Rq=8 kpc and 
Zq=15 pc l)Yamagatnr^bshiil 119921) are the solar ra- 
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Fig. 3. — A sample realization of the Monte Carlo simulation with 
10 10 random points with the best-fit vertical scale height of 350 pc. 
The over-plotted crosses with numbers represent e ach of the ob- 
served HST/ACS field in Table |T] The effects of the lSchlegel et alJ 
( 1998) dust maps are readily apparent when comparing the upper 
(no extinction) and lower (extinction) panels. 

dius and height, respectively. The vertical scale height, 
Zh, is the free parameter and found by minimizing the 
squared difference between the number counts from the 
model and the HST data. Altering the assumed coordi- 
nates of the Sun and the radial scale length have little 
effect on the vertical scale height estimate. To generate 
absolute m agnitudes, we adop ted the J-band luminosity 
function oflCruz et aU l)2003|) and the (V — J) colors of 
IHawlev et alJ l|2002|) over the appropriate range of spec- 
tral type. 

The effects of interstellar extinction wer e included into 
the m odel using the COBE dust maps of iSchleeel et al.l 
(1998) in two ways. First, we assumed that each point 
was located beyond the dust, establishing a lower bound 
on the model counts. Alternatively, an upper bound is 
reached by assuming that the Galaxy has no dust what- 
soever. Since an overwhelming majority of the dust is 
localized to \b ll \ < 15° (see Figure|3J) and only two of our 
observed ACS fields are in this range, either approach 
yielded the same result within the uncertainties, there- 
fore we adopted the third method for simplicity. Fig- 
ure |3] is a representative realization the model with a 
scale height of 350 pc, the with locations of the 15 ob- 
served fields are indicated with plus signs. 

The canonical di sk/spheroid G alaxy likely has addi- 
tional components (Bahcall 1986), the models used here 
did not contain any contr ibution from the Galactic bulge 
or a t wo-component disk ( Gilmqre & Reid 1983i lGilmorei 
1984) for the following reasons. First, a bulge distribu- 
tion was not modeled since its radius is ~l-2 kpc or 
~7°-14°, and every field is well beyond 14° of the Galac- 
tic center, hence we do not require a bulge component 
in the model. Second, this sample contains only 28 L- 
& T-dwarfs which are likely within ^1000 pc (based on 
the luminosity function) of the Sun. Since the thick disk 
has a scale height of >1000 pc, we expect the star counts 
to be dominated by a single disk population. Moreover, 
With only 28 candidates the models and analyses must 
remain simple and straightforward. 

4. ANALYSIS 

Despite this work utilizing the largest dataset of L- 
& T-dwarfs compiled from HST observations, the star 
counts remain very sparse requiring a simple analysis 
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Fig. 4. — Top panel: model surface density as a function of 
Galactic latitude. Here we have averaged over all longitudes with 
a 2.5a--clipping to better handle the sparse statistics at high lati- 
tudes. The dotted and full lines indicate models with and without 
the extinction corrections, respectively. The data from the 15 ACS 
fields are plotted for comparison as filled circles and downward ar- 
rows as an upper limit when either zero or one object was detected. 
Bottom panel: the residuals from the upper panel as a function 
of Galactic latitude. Clearly the data where the dust corrections 
are large (ie. |b n | <15°) are the most deviant. The model used in 
both panels has a vertical scale height of 35 pc . The two fields 
outlying fields at \b u \ <15° are discussed in § 14.11 

scheme. Using the grid of Monte Carlo models described 
in § [3J we sought the scale height which minimizes the 
squared difference between the integrated star counts of 
the model and those from the HST/ACS dataset. For 
the fields where no L- & T-dwarf candidates were found, 
we assumed an upper limit of one object (per field) could 
have been detected and perform this minimization tech- 
nique simultaneously on all 15 fields. This procedure 
yielded a vertical scale height of 350±50 pc. In the upper 
panel of Figure^] we plot the modeled surface density av- 
eraged over Galactic longitude as a function of Galactic 
latitude as computed from the model with a scale height 
of 350 pc, with the HST data points from Table over- 
plotted for comparison. The residuals in the lower panel 
clearly demonstrate that the model with a scale height 
of 350 pc reproduces the HST star counts for |fr n | > 15° 
where dust extinction is minimal. 

4.1. Sources of Contamination 

While all L- & T-dwarf candidates were systemati- 
cally found by color and FWHM criteria, each was vi- 
sually confirmed as a point source. However, the color- 
magnitude rules outlined in § potentially find three 
classes of contaminates: 

(1) The primary motivation of this study is to re- 
liably correct the z ~ 6, i'-band dropo ut galaxy 
surveys for interloping L- fc T-dwarfs iYan et al 



l2002llYan. Windhorst, fc CoheiJ200.3IYan fc Windhorst 



Since our method is similar to the i'-band 
dropout technique, we expect possible contamination 
to the star counts from the z ~ 6 galaxies. In a re- 
cent study of the Hubbl e Ultra-Deep Field (HUDF), 
lYan fc Windhorst! l|2004bD find 108 i'-band dropouts us- 
ing the {%' — z')>1.3 mag color selection. Their sam- 
ple has a median z'-band magnitude of 28.5 mag, and 
contains only three objects brighter than our limit of 
z'=26.0 mag. Each of these three objects is considerably 
extended and could not be mistaken for a point-source. 
Since the majority of unresolved z ~ 6 galaxies will be 
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~2.5 mag too faint, we conclude their contamination in 
our sample must be negligible. 

(2) Another known source of possible contaminates 
comes from dusty, elliptical galaxies with redshifts 1.0 < 
z < 1.5, whose 4000 A break occurs betwe en the i' and z' 
bands fYan. Windhorst, fc Cohenl I2003H . With a typ- 
ical color of (i' — z')~1.0 mag, many of these objects 
would appear too blue in the absence of extreme internal 
reddening. Moreover the visual identification confirms 
only point-like objects populate the "stellar locus" in the 
lower left of Figure [21 Hence, extended elliptical galaxies 
could not grossly corrupt our sample. 

(3) In addition to the above extragalactic sources, we 
anticipated contamination from galactic M-dwarfs for 
two different reasons. First, the color criterion of (i'—z')> 
1.3 mag was primarily motivated by the z ~ 6 galaxy s ur- 
veys and is ~0.5 mag too blue (Ha wlev et al.1 12002]) to 
have included only L- & T-dwarfs. When we repeated the 
above analysis for (i'— z')>1.8 mag (Ha wlev et al.11200'2^ 
the star counts were reduced by ~50% and the inferred 
vertical scale height was 300 ±100 pc. While there is a 
significant contribution from late M-dwarfs, the vertical 
scale height was unaffected by the 0.5 mag color differ- 
ence. Second, an appreciably reddened M-dwarf could 
have an (i' — z') color of an unreddened L- or T-dwarf. 
This scenario requires considerable reddening, only the 
two fields with |6 n |<15° have E(i'-z') >0.1 mag. This 
effect could account these fields lying more than la above 
the best-fit line in Figure^ however we cannot be certain 
without more broad-band filters or spectroscopy. We in- 
vestigated this effect's contribution by removing the two 
suspect fields and repeated the analysis. While this pro- 
cedure reduced the observed star counts by ~50%, it 
resulted in a vertical scale height of 360±180 pc. With- 
out further observations it is difficult to definitively re- 
move highly reddened M-dwarfs, however their contribu- 
tion should not grossly affect our main goal. 



5. DISCUSSION 

Using a suite of Monte Carlo simulations and 15 
HST/ACS parallel fields, we find a vertical scale height 
of 350±50 pc for the L- & T-dwarf population based on 
28 faint candidates. This estimated scale height is consis- 
tent with the known trend of increasing scale height with 
decreasing stellar mass, independent of reddening, color 
selections, and other Galactic parameters and is within 
the uncertai nties of and is l argely a refinement of pre- 
vious work ijLiu et al.l l2002: Pir zkal et al]l2005|) . Using 
our value of the scale height, and the parameters given 
in §|21 we predict a total of ~10 n L- & T-dwarfs and a 
total mass of < 10 9 M Q in the Milky Way. 

This improved understanding of the L- & T-dwarf 
Galactic distribution will aide high-redshift surveys 
in better estimating the contamination of L- & T- 
dwarfs in their samples. In the recent HUDF 
pointed at (a, <5)=(3 h 32 m 39!0,-27°47'29('l) and a depth 
of z 1 ~29 mag, we predict >2 L- & T-dwarfs in its 
~11 arcmin 2 field o f view, which has been confirmed by 
iPirzkal et alJ l)2005l) who have spectroscopically identi- 
fied three L- & T-dwarfs. We confirm that Galactic L- 
& T-dwarfs cannot significantly corrupt the z ~ 6 surveys 
in high-latitude fields (the HUDF for example), however 
low-latitude fields will find a modest number of inter- 
loping L- & T-dwarfs. With only 28 candidates from 
15 fields, our statistics remain sparse and ideally require 
further observations. 
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